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GLACIATION IN ELYSIUM. Duwayne M. Anderson, 
Texas A&M University, College Station TX 77843-3115, 
USA. 

Results of a study of high-resolution Viking Orbiter 
images of the northwestern slopes of Elysium Mons, utilizing 
a variety of image analysis techniques, provide striking 
evidence of an extended period of glaciation that involved a 
large ice sheet of greater than 1 .2 km thickness and subice 
fluvial activity that contributed to the formation of the 
channels of Hrad Valles and Granicus Valles. These two 
unusual channel systems begin on the lower slopes of 
Elysium Mons and extend into Utopia Planitia in the region 
lying between 215° and 230°W and between 25° and 45°N. 
Key indicators are the presence of serrated volcanic 
constructs whose dimensions and morphology indicate an 
origin involving successive, localized, subice fissure 
eruptions. The channels visible in Hrad Valles bear a striking 
resemblance to subice fluvial features found in the dry valleys 
of Antarctica at the periphery of this massive continental ice 
sheet. Crater size distributions and crater morphologies are 
quite consistent with this interpretation, as are certain other 
topographic features suggesting the presence in the past of 
ice-rich permafrost that, while having undergone significant 
degradation, still remains. 
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ORBITAL, ROTATIONAL, AND CLIMATIC INTER- 
ACTIONS: LESSONS FROM EARTH AND MARS, 

Bruce G. Bills, Geodynamics Branch, NASA Goddard Space 
Flight Center, Greenbelt MD 20771, USA. 

Introduction: Though variations in orbital and rota- 

tional parameters of the Earth and Mars are widely rec- 
ognized as plausible sources of significant climatic variation 
on 10 3 -10 8 -yr timescales, many aspects of the connection 
between orbital, rotational, and climatic variations remain 
poorly understood. In general, the orbital histories are very 
well known, the rotational histories are less well known, and 
the climatic histories (especially for Mars) are very poorly 
known. A brief review will be given of recent progress in 
computing orbital and rotational secular variations, and in 
connecting them to climatic change. The emphasis will be on 
highlighting those areas that limit our present understanding. 

It is obvious that mass redistributions associated with cli- 
matic change (glaciations) are an important source of crustal 
deformation and geodynamic change on the Earth, and may 
have played a similar role on Mars in the distant past. It is 
much less widely appreciated, however, that rates, phases, 
and amplitudes of deformation of the deep interior of the 
planet can influence climate. The mantle and core, if com- 
pletely decoupled, would precess at quite different rates, and 
even with plausible coupling strengths, some degree of dif- 
ferential precession is possible. Mass flow associated with 


glaciations can influence the processional dynamics of the 
Earth or Mars, and changes in orbital and rotational parame- 
ters influence the latitudinal and seasonal pattern of insola- 
tion. Previous attempts to account for astronomically forced 
climatic change have usually only considered extremely 
simplisitic models for the response of the solid planet to ex- 
ternal torques and surface loads. 

The Traditional Perspective: The latitudinal and sea- 

sonal pattern of incident solar radiation depends on the ec- 
centricity of the planetary orbit and the orientation of the spin 
axis relative to both the orbit normal and the aspidal line. 
Unit vectors s and n characterize the directions of the spin 
axis and orbit normal respectively. Two angles completely 
characterize the relative orientation of the spin axis. The 
obliquity e is simply the angle between the orbit normal and 
the spin axis 


e = cos 1 (n*s) (1) 

The ascending node of the orbit plane on the instantane- 
ous equator plane has an orientation given by (s x n), and the 
longitude of perihelion to is just the angle in the orbit plane 
from that node to perihelion. It is widely appreciated that 
secular variations in these three parameters (e,e,to) produce 
major climatic change [1,2]. In fact, spectral analyses of long, 
high-resolution marine sediment isotopic records show sig- 
nificant variance at periods near 100 kyr, 41 kyr, and 
19-23 kyr, which are generally attributed to spectral lines in 
the radiative forcing fluctuations associated with e, e, and e 
sin( to) respectively. 

The causes and effects of the orbital changes are quite 
well understood. Gravitational interactions with the other 
planets cause the shape and orientation of the orbit to change 
on timescales of 10 4 — 10 6 yr. The inclination I and nodal lon- 
gitude Q determine the orientation of the orbit plane. The ec- 
centricity e and perihelic longitude co determine the shape of 
the orbit and its orientation within the plane. Note that © is 
measured from an inertially fixed direction, rather than the 
moving node as is the case for tu. The secular evolution of the 
orbital element pairs (1,0) and (e,co) can be conveniently rep- 
resented in terms of Poisson series 


p = sin(I)sin(0) = ZNjSin(sj t + gj) 
q = sin(I)cos(Q) = ENjcos(sj 1 4* gj) 


h = esin(co ) = IMj sin(rj t + fj) 
k = ecos(co) = EMj cos(rj t + fj) 


( 2 ) 

( 3 ) 


In the lo west-order solution, there are as many frequencies 
r and Sj as there are planets. However, the frequencies and 
Sj are characteristic modal frequencies (eigenvalues) of the 
coupled system of oscillators and are not each uniquely asso- 
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dated with a particular planet [3], The frequencies are all 
positive, indicating that the perihelia advance. In the lowest- 
order solution, the apsidal rates are all in the interval 
(0,667 < t- < 28.221 arcsec/yr). The corresponding periods 
are 45.92 kyr to 1.943 m.y. One of the frequencies, s jf is zero, 
and all the others are negative, indicating that the nodes re- 
gress. In the lowest-order solution, the nonzero nodal rates 
are all in the interval (0.692 < Sj < 26.330 arcsec/yr). The 
corresponding periods are 49,22 kyr to 1.873 m.y. In higher- 
order solutions, variations in (e,co) become coupled to vari- 
ations in (I,Q), but the solutions can still be cast in terms of 
Poisson series like equations (2) and (3). 

Laskar [4] has recently published a secular variation 
theory that is complete to fifth order in eccentricity and 
inclination. Agreement between this secular variation model 
and strictly numerical computations [5,6] is much better than 
for any previous analytical model. The inclination and 
eccentricity series for Earth and Mars each contain 80 distinct 
terms. 

In computing these secular orbital variations, the Earth, 
Moon, and planets can all be treated as point masses. No in- 
ternal structure or processes are relevant to orbital evolution. 
The physics of the process is simple and well understood, 
though development of proper mathematical tools to repre- 
sent the long-term evolution remains an area of active re- 
search [4,6,7]. On the other hand, the rotational evolution 
does depend rather sensitively on various aspects of the struc- 
ture and dynamics of the interior. 

Lunar and solar gravitational torques acting on the oblate 
figure of the Earth cause the spin axis s to precess about the 
instantaneous orbit normal n. If the Earth is considered to be 
a rigid body, the evolution of the spin axis orientation is 
given by 

dsdt = a(n • s)(s x n) (4) 

where 

(5) 

is a scalar rate factor that depends on intrinsic properties of 
the Earth, such as polar and equatorial moments of inertia 
(CA) and rotation rate n, and on extrinsic influences, such as 
masses m, orbital inclinations I, and semiminor axes b, of the 
Moon and Sun. The solar and lunar torques together produce 
a precession of the spin axis of the Earth at a rate of a(n • 
s) = 50.38 arcsec/yr [8,9]. 

Unfortunately, in the case of Mars, the precession rate is 
still rather poorly known. The best estimates at present come 
from analysis of the Viking Lander range data, which yield 
values of 9.6 ± 0.6 arcsec/yr. Much smaller relative errors are 
frequently cited, but all such optimistic estimates are directly 


dependent upon an assumed value for the moment of inertia 
of Mars, which is highly model dependent [10], precisely be- 
cause the axial precession rate is not known. The uncertainty 
in this parameter is the largest single impediment to accurate 
reconstructions of the obliquity history of Mars [11]. Fortu- 
nately, the Mars Observer mission radio science investiga- 
tions will significantly improve our knowledge of the pre- 
cession rate of Mars within a few months after orbit inser- 
tion. 

Once the present spin axis direction s is known and 
orbital element histories are given via equations (2) and (3), 
an obliquity history can be constructed from equation (4) in 
two different ways. The linear perturbation approach [12-17] 
involves deriving coefficients of a trigonometric series, simi- 
lar to equations (2) and (3), which yields the obliquity and 
longitude of perihelion directly as functions of time. An al- 
ternative is to apply standard numerical algorithms for solv- 
ing initial value problems to generate a vector time series s(t) 
and then compute the obliquity and longitude of perihelion 
directly [11,18,19], The spectrum of obliquity variations, in 
the linear perturbation model, is simply obtained from the 
inclination spectrum by shifting each frequency s } by the 
1 uni solar precession rate (a = 50.38 arcsec/yr for the Earth, or 
8-10 arcsec/yr for Mars) and multiplying each amplitude Nj 
by the spectral admittance 



There are several important features of this solution: 
(l)As a is positive, and s } is negative in all cases, it is 
possible for the denominator to vanish; (2) near this 
resonance, the obliquity variations are large; and (3) at 
smaller or larger forcing frequencies, the obliquity variations 
are small. 

Uncertainty in the obliquity history of Mars derives from 
two facts: the value of the precession rate is uncertain by 6% 
(compared to 0.04% for the Earth) and the obliquity history 
of Mars depends more sensitively on the precession rate than 
is the case for Earth because of a near resonance with some of 
the inclination forcing frequencies. To investigate the sensi- 
tivity of the computed martian obliquity history to assumed 
precession rate, a series of numerical integrations of the rigid 
body precession equations was made [11], covering a 20-m.y. 
interval centered on the present, using precession rates in the 
interval 8. 5-9.5 arcsec/yr. The maximum obliquity encoun- 
tered can be anywhere in the range 35°-50°, and the 
minimum value can be 8M4®. Despite this large uncertainty 
in the particulars of the obliquity history, all the computed 
time series were characterized by two dominant features: 
large oscillation with a characteristic period of ~10^ yr and a 
significant modulation with a characteristic period of 2 
10 6 yr. 
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Fourier, Legendre, and Milankovitch: The influence 

of orbital and rotational variations on climate is operative 
through perturbations in the latitudinal and seasonal pattern 
of insolation. The diurnal average intensity of radiation at a 
point is inversely proportional to the squared solar distance 
and directly proportional to the diurnal average rectified solar 
direction cosine 

F = (a/r) 2 (||u-u,||) (7) 

where a and r are mean and instantaneous solar distance, and 
u and u s are unit vectors from the center of the Earth to the 
surface point of interest and the subsolar point respectively. 
The insolation pattern, as a function of latitude 0 and mean 
anomaly M, can be readily computed once values are speci- 
fied for the orbital and rotational parameters c, e, and w 
[14,15,18,20,21]. This pattern can also be written in terms of 
a Fourier-Legendre series [20,22-24] 

F0i,M;e,e,ro) = £P„(n )£exp(ipM) T n Je,e,w) (8) 

where Ji = cos(0) and P n is a Legendre polynomial. The 
number of terms in the Fourier summation required to obtain 
a good representation of the seasonal pattern is greater in the 
polar regions than in the tropics and mid latitudes. The pri- 
mary difficulty in the polar regions is reproducing the abrupt 
change in slope of the insolation curve at times of transition 
to continual darkness or continual light. It is also true that the 
polar regions place the greatest demands on the Legendre 
summation, since the spatial pattern also has a discontinuous 
first derivative at the latitude where the transition occurs to 
continual darkness or light. 

Processional Dynamics with Variable Rate: All but 

the most recent reconstructions of the radiative forcing input 
to paleoclimate models have assumed that both the orbital 
and rotational dynamics could be readily and accurately 
reconstructed from their present configurations via the simple 
analyses mentioned in the introduction. These expectations 
seem well founded in the case of orbital evolution, though the 
possibility of chaotic dynamics in the inner solar system 
[7,25] does seem to preclude confident extrapolation beyond 
10 7 yr. However, there are a number of processes, working in 
different locations and at different rates, that all serve to 
compound the difficulty of accurately computing the spin pro- 
cessional evolution . 

On the longest timescales of interest (10 7 ~10 9 yr) the 
limiting uncertainty is variability in the tidal transfer of 
angular momentum from the rotation of the Earth to the orbit 
of the Moon. At present, these tidal torques are increasing the 
length of the day by 22.5 1CT 6 s/yr and increasing the size of 
the lunar orbit by 3.88 cm/yr [26,27]. Berger et al. [28] have 
made a useful first step toward including this effect in cli- 


matic time series. They computed the change in the major 
precession and obliquity frequencies due to lunar tidal evolu- 
tion assuming that the present rate of tidal energy dissipation 
is representative of the past 500 m.y. However, the present 
rates are considerably higher than the long-term average [29], 
largely due to a near resonance between sloshing modes of 
ocean basins and the diurnal and semidiurnal tidal periods 
[30], and apparently compounded by a contribution from 
shallow seas [31,32]. Sedimentary records that constrain lu- 
nar orbital evolution show some promise of resolving this 
problem [33-36], but the situation is definitely more complex 
than is suggested by Berger et al. [28]. 

Another parameter that can vary, on rather shorter time- 
scales and in an equally irregular fashion, is the gravitational 
oblateness of the Earth (C-A)/C. Thomson [37] has recently 
made three important contributions to the understanding of 
this source of variability. First, he pointed out that mass re- 
distribution associated with major glaciations and compensat- 
ing subsidence and crustal deformations [38,39] can cause 
fractional changes in oblateness of order 10 -3 -10 -2 . Second, 
he showed that high-resolution spectral analyses of several 
climatic time series appear to indicate fluctuations of the 
luni solar precession rate of this magnitude, and with a 
dominant period near 100 kyr. Finally, Thomson pointed out 
that the best fit to the paleoclimate proxy data was obtained 
using a mean luni solar precession rate 0.6 arcsec/yr less than 
the present observed value. He notes that the resulting value 
would correspond rather closely with that expected for a hy- 
drostatic flattening [40]. If these important results are cor- 
roborated, they will demonstrate that important feedback 
loops exist in the orbital-rotational-climatic interactions sys- 
tem, further “up-stream” in the presumed causal chain than 
has been previously recognized. 

Differential Precession of the Mantle and Core: The 

hydrostatic figure of a planet represents a compromise be- 
tween gravitation, which attempts to attain spherical symme- 
try, and rotation, which prefers cylindrical symmetry. Due to 
its higher mean density, the core of the Earth is more nearly 
spherical than the mantle. The direct lunisolar processional 
torques on the core will thus be inadequate to make it process 
at the same rate as the mantle. In fact, the core oblateness is 
only about three-fourths that required for coprecession with 
the mantle [41]. However, it is clearly the case that the core 
and mantle precess at very nearly the same rate [42]. A va- 
riety of different physical mechanisms contributes to the 
torques that achieve this coupling, but a purely phe- 
nomenological partitioning is useful. The net torque can be 
described as a sum of inertial torques, which are parallel to 
(Xm x Xc)> and dissipative torques, which are parallel to ~ 
Xc)- Here, Xc and Xm are the rotation vectors of the core and 
mantle respectively. The two types of torques have qualita- 
tively different results: Inertial torques cause the core and 
mantle axes to precess at fixed angular separations and on the 
opposite side of their combined angular momentum vector, 
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whereas the effect of dissipative torques is to reduce the 
angle between the axes. 

On short timescales it is appropriate to consider the core 
to be an in viscid fluid constrained to move within the 
ellipsoidal region bounded by the rigid mantle [43-45]. The 
inertial coupling provided by this mechanism is effective 
whenever the ellipticity of the container exceeds the ratio of 
the precessional to rotational rates. If the mantle were 
actually rigid, or even elastic [46,47], this would be an 
extremely effective type of coupling. However, on sufficiently 
long timescales, the mantle will deform viscously and can 
accommodate the motions of the core fluid [48]. The inertial 
coupling torque exerted by the core on the mantle will have 
the form 

Tj = ki[X m x X c ] (9) 

A fundamentally different type of coupling is provided by 
electromagnetic or viscous torques [49-51]. The dissipative 
coupling torque exerted by the core on the mantle will have 
the form 

Ta=kd[XmxXc] (10) 

This type of coupling is likely to be most important on 
longer timescales. In each case, the mantle exerts an equal 
and opposite torque on the core. The response of the coupled 
core-mantle system to orbital forcing is given by [52-54] 

ds^/dt = O m (n- sJ(Sm *n)-Pm(Sm- s c ) - Ym( s m * *) 
dSc /dt = o^(n ■ s^X^c x n) + PcCs,,, - s c ) + y c (s m x s c ) 

where a m is similar to a above, except that only mantle mo- 
ments A m and C m are included, and 

( 12 ) 

Ym=ki/C m o 2 

where o is the mean rotation rate. 
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CONDENSATION PHASE OF THE MARTIAN SOUTH 
POLAR CAP. J. Capuano, M. Reed, and P. B. James, De- 
partment of Physics and Astronomy, University of Toledo, 
Toledo OH 43606, USA. 

One type of database that can be useful in constraining 
models of the martian surface-atmosphere system is the time- 
dependent boundary of CO, frost for the polar caps. These 



